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ABSTRACT 

The study of the magnetic field in stellar radiation zones is an important topic in 
modern astrophysics because the magnetic field can play an important role in several 
transport phenomena such as mixing and angular momentum transport. We consider 
the influence of rotation on stability of a predominantly toroidal magnetic field in 
the radiation zone. We find that the effect of rotation on the stability depends on 
the magnetic configuration of the basic state. If the toroidal field increases sufficiently 
rapidly with the spherical radius, the instability cannot be suppressed entirely even 
by a very fast rotation although the strength of the instability can be significantly 
reduced. On the other hand, if the field increases slowly enough with the radius or 
decreases, the instability has a threshold and can be completely suppressed in rapidly 
rotating stars. We find that in the regions where the instability is entirely suppressed a 
particular type of magnetohydrodynamic waves may exist which are marginally stable. 

Key words: instabilities - magnetohydrodynamics - stars: interiors - stars: magnetic 
field - Sun: interior 



1 INTRODUCTION 

In recent years, magnetic fields of various strength and 
topology have been detected in increasing number of stars 
(see, e.g., Donati et al. 2006). Despite the origin of stellar 
magnetic fields is usually related to the turbulent dynamo 
action in convection zones, it is quite possible that magnetic 
fields exist also in internal radiation zones as well. For in- 
stance, a uniform rotation of the solar core may be explained 
naturally by the presence of the magnetic field. A thin transi- 
tion layer between the convection and radiation zones of the 
Sun, called the tachocline, cannot be understood within the 
frame of purely fluid-dynamical mechanisms and requires a 
large-scale magnetic field in the Sun's interior (Gough & 
Mclntyre 1998). The origin, topology, and strength of the 
magnetic field in internal radiation zones are the subject of 
debates for decades. Likely, dynamo cannot operate in stellar 
radiation zones. A dynamo action requires sufficiently rapid 
flows with Re m 3> I (Re m is the magnetic Reynolds num- 
ber) to generate a global magnetic field but, according to the 
generally accepted point of view, such flows are not avail- 
able in the internal radiation zones (see, e.g., Schwarzschild 
1958). Sometimes, the origin of the magnetic field in radia- 
tion zones is related to the dynamo mechanism proposed by 
Spruit (1999) but the existence of such mechanism has not 
be proven yet. One might speculate that a fossil magnetic 



field could exist in radiation zones. Even a weak fossil field 
with non-vanishing poloidal component will quickly wrap 
up into a predominantly toroidal configuration, under the 
action of differential rotation. Such configurations can be 
generated if Re m of differential rotation is greater than 1, or 
[Vfi| > r/ m /r 3 where ft and r/m are the angular velocity and 
magnetic diffusivity, respectively. Estimating |Vfl| ~ AQ/r 
where Afi is a departure from the rigid rotation and assum- 
ing that the conductivity of plasma is ~ 10 16 s _1 , one can 
obtain that this condition is satisfied if Afi/n > lO -18 ^^ 
where n sec is the angular velocity in inverse seconds. There- 
fore, even very weak departures from the rigid rotation lead 
to a generation of the toroidal field. The magnetic field 
with a predominant toroidal component is also typical for 
the liquid cores of neutron stars (Bonanno et al. 2005, 2006) 
where the field can be generated by the turbulent dynamo 
during the very early evolutionary stage. The toroidal field 
resulting from winding up or dynamo action cannot be arbi- 
trarily strong because various magnetic instabilities will set 
in when it becomes strong enough. That is why much atten- 
tion has been drawn to the instabilities that may affect the 
toroidal magnetic field in stably stratified radiation zones. 

The magnetic field in a radiation zone can be subject 
to various instabilities such as the magnetic buoyancy (see, 
e.g., Oilman 1970, Acheson 1978) or magnetorotational in- 
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stability (Velikhov 1959, Balbus 1995). Likely, however, that 
the most efficient instabilities are caused by the electric cur- 
rents maintaining the magnetic configuration (Spruit 1999). 
Such instabilities are well studied in cylindrical geometry in 
the context of laboratory fusion research (see, e.g., Freid- 
berg 1970, Goedbloed & Hagebeuk 1972). In astrophysical 
conditions, the instability caused by electric currents is stud- 
ied mainly in cylindrical geometry as well (see, e.g., Tayler 
1973). It turns out that the properties of instability depend 
on the ratio of the axial and toroidal fields and, even if this 
ratio is small, the axial field can alter the instability substan- 
tially (see, e.g., Knobloch 1992, Bonanno & Urpin 2008a). 
The effect of an axial field on the Tayler instability of the 
toroidal field has been studied in detail in cylindrical geom- 
etry by Bonanno & Urpin (2008b, 2011). The nonlinear evo- 
lution of the Tayler instability was considered by Bonanno et 
al. (2012) who argued that symmetry-breaking can give rise 
to a saturated state with non-zero helicity even if the initial 
state has zero helicity. A production of non-zero helicity is 
likely possibility for the onset of adynamo action. 

Stability of the spherical magnetic configurations is 
much less studied because of mathematical problems. With 
numerical simulations Braithwaite & Nordlund (2006) stud- 
ied the stability of a random initial field and argued that 
this field relaxes on a stable mixed magnetic configuration 
with both poloidal and toroidal components. Stability of 
magnetic configurations with a predominantly toroidal field 
draws a particular attention. Numerical modeling by Braith- 
waite (2006) confirmed that the toroidal field with B v oc s 
or oc s 2 (s is the cylindrical radius) is unstable to the m = 1 
mode as it was predicted by Tayler (1973). It is widely be- 
lieved that rotation and stratification can play an important 
role in stellar radiation zones providing a stabilizing influ- 
ence on the Tayler instability (see, e.g., Pitts & Tayler 1985, 
Spruit 1999). Recently, the effect of stratification and ther- 
mal conductivity on the Tayler instability of the toroidal 
field has been considered by Bonanno & Urpin (2012). The 
authors studied a linear stability and used a local approxi- 
mation in latitude and global in the radial direction. They 
calculated the growth rate of instability and argued that 
the combined influence of gravity and thermal conductiv- 
ity can never suppress the instability entirely. Stratification 
suppresses the instability at the pole more efficiently than at 
the equator. The growth rate of instability can be essentially 
reduced by a stable stratification. A decrease of the growth 
rate caused by stratification is inversely proportional to the 
Brunt- Vaisala, frequency. Therefore, if gravity is strong, the 
instability developes very slowly. A simple fitting expression 
has been obtained for the growth rate of instability in a 
stratified radiation zone. Bonanno & Urpin (2012) argued 
that the instability of modes with a large number of nodes 
in the radial direction is significantly more suppressed than 
the fundamental eigenmode. The reason of this is qualita- 
tively clear. As it was pointed out by Spruit (1999), the 
stabilizing influence of stratification is less pronounced for 
perturbations with short radial lengthscales. Therefore, it 
seems at the first glance that the instability should operate 
most efficiently on very short radial lengthscales. However, 
this conclusion is incorrect because a destabilizing effect of 
electric currents (that is the reason of the Tayler instability) 
also decreases if the radial lengthscale decreases. This occurs 
because the instability becomes almost two-dimensional if 



the radial lengthscale goes to but, as it was shown by 
Kitchatinov & Riidiger (2008), the Tayler instability does 
not exist in 2D. Therefore, the instability cannot arise for 
perturbations with very short radial lengthscales despite the 
stabilizing influence of stratification is minimal for them. 

Rotation can also suppress the Tayler instability and 
stabilize the magnetic configurations. The effect of rotation 
has been considered by a number of authors. Spruit (1999) 
found that the growth rate of the Tayler instability near 
the rotation axis should be of the order of ~ ua(wa/!!) if 
£2 ^> u>a in a particular case of the magnetic field depen- 
dent on s alone, where lo a and Q are the Alfven frequency 
and angular velocity of the star, respectively. In numerical 
modelling by Braithwaite (2006), the Tayler instability was 
suppressed if Q is above a certain value of the order of uja- 
Above this value a distinct oscillatory behaviour sets in with 
marginal stability. Stability of the toroidal field in rotat- 
ing radiation zones has been considered also by Kitchatinov 
(2008) and Kitchatinov & Riidiger (2008) who clarified that 
the Tayler's instability recovers only in 3D. Radial displace- 
ments are essential for this instability. It does not exist in 
the 2D case of strictly horizontal (perpendicular to gravity) 
disturbances and only stable modes exist in this case. These 
authors argued that the magnetic instability is character- 
ized by the threshold field strength. The instability arises if 
the magnetic field is stronger than this threshold but it does 
not occur if the field is weaker than the threshold. Estimat- 
ing this threshold in the solar radiation zone, Kitchatinov 
& Riidiger (2008) impose the upper limit on the magnetic 
field « 600 G. The stability of the toroidal field in a rotating 
radiation zone has been studied also by Zahn et al. (2007) in 
the particular case B v oc s. The particular type of oscillatory 
modes found by these authors is relevant to rotation and is 
stable in the non-dissipative limit. However, instability of 
the considered modes can occur in a form of an oscillatory 
diffusive instability if dissipation is provided by radiative or 
Ohmic diffusion. Unfortunately, the authors of these stud- 
ies did not compare the results of their calculations and the 
reason of such qualitatively different behaviours is unclear. 

In this paper, we consider in detail the effect of rota- 
tional suppression of the Tayler instability in the case of 
a predominantly toroidal field. We show that rotation can 
influence the instability in different ways depending on the 
baseground magnetic configuration. The paper is organized 
as follows. The basic equations and mathematical formula- 
tion of the problem are presented in Sec. 2. This is followed 
by results of numerical calculations of the growth rate and 
frequency of unstable modes in Sec. 3. The paper closes with 
a summary of the main results and some remarks in Sec. 4. 



2 BASIC EQUATIONS 

We assume that the magnetic field in a radiation zone is sub- 
thermal and the magnetic pressure is smaller than the gas 
pressure. In this case, the Boussinesq approximation is ap- 
plied for a consideration of low-frequency modes. The ideal 
MHD equations read in this approximation 

^ + (a-\7)v = -^+g+-L(X7xB)xB, (1) 
ot p 4irp 
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V-w = 0, V-B = 0, 



(2) 
(3) 



where g is gravity, p and p are the gas pressure and density, 
respectively. The equation of thermal balance reads in the 
Boussinesq approximation 



^+v-(VT-V ai T) = Q, 



(4) 



where V a <jT is the adiabatic temperature gradient. 

Consider the stability of an axisymmetric toroidal mag- 
netic field using spherical coordinates (r, 9, tp). We assume 
that the radiation zone rotates with the angular velocity 
0=const and that the toroidal field depends on r and 9, 
B v — B v (r,9). In the unperturbed state, the radiation zone 
is assumed to be in hydrostatic equilibrium, then 



— = g + -J—CV x B) x B + e s Q, 2 
p 4-zrp 



(5) 



where e s is the unit vector in the cylindrical radial direc- 
tion. The rotational energy is assumed to be much smaller 
than the gravitational one, g ^> rQ. 2 . Since the magnetic 
energy is subthermal, g is approximately radial in our basic 
state. Therefore, longitudinal variations of the unperturbed 
density and pressure are small. 

We consider a linear stability. Small perturbations will 
be indicated by subscript 1, while unperturbed quantities 
will have no subscript. Linearizing Eqs.(l)-(4), we take into 
account that small perturbations of the density and tem- 
perature in the Boussinesq approximation are related by 
pi/p = —fi(T\/T) where /3 is the thermal expansion co- 
efficient. We use a local approximation in the (9-direction 
and assume that small perturbations depend on 9 as oc 
exp(— U9), where I ^> 1 is the polar wavenumber. Since the 
basic state is stationary and axisymmetric, the dependence 
of perturbations on t and tp can be taken in the exponen- 
tial form as well. Then, perturbations are proportional to 
exp (at — U9 — imtp) where m is the azimuthal wavenumber 
and a is the growth rate. The dependence of perturbations 
on r should be determined from Eqs.(l)-(4). 

Since the influence of gravity on instability has been 
studied in detail by Bonanno & Urpin (2012), we concentrate 
in this paper mainly on the effect of rotation. We consider 
a simplified problem assuming that stratification is neutral 
and VT = X7 a dT. In this case, the stabilizing effect of strat- 
ification is excluded and we can study only the influence 
of rotation alone. The combined influence of rotation and 
stratification will be considered elsewhere. For the sake of 
simplicity, we also assume that the unperturbed density is 
approximately homogeneous in the radiation zone. 

Eliminating all variables in favor of vi r , we obtain with 
the accuracy in terms of the lowest order in (kgr)^ 1 



^2 2 2 2 22^2/^1 1 

r 2 r \H k 2 

2 4° 



(6) 



laoile 1- — =- 

r rk v <7q 



vi r = 0, 



where the prime denotes a derivative with respect to r and 

D = . (7) 
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a — vmil, lj a - 
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tti = 20, cos 9, Q e = 2Q sin 9, 



0% +UJ 
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k± — k$ + ku^, 



H = 9- r HrB v ). 



The polar and azimuthal wavevectors are kg = l/r and k v — 
rajr sin 9, respectively. 

This equation with the corresponding boundary condi- 
tions describes the stability problem as a non-linear eigen- 
value problem. Fortunately, the main qualitative features 
of this problem are not sensitive to the choice of boundary 
conditions. That is why we choose the simplest boundary 
conditions and assume that v\ r = at the inner and outer 
boundaries, r = Ri and r — R, respectively. 



3 NUMERICAL RESULTS 

We assume that the radiation zone is located at Ri r ^ R 
or, introducing the dimensionless radius x = r/R, at Xi ^ 
x ^ 1 where Xi = Ri/R. We choose the internal radius of the 
radiation zone, Xi, to be small but finite for computational 
convenience. In particular Xi = 0.1 in this investigation. Our 
results do not depend qualitatively on the precise value of 
Xi, as we have explicitly verified by performing calculations 
with with Xi in the range 0.1 < Xi < 0.4 (less than 1% 
difference in the growth rate for Xi = 0.1 and Xi — 0.3 for 
a = 2 for instance) . 

We represent the toroidal field as 



B v = Bi)(x/xi) a sin#, 



(8) 



where Bo is the field strength at x = Xi at the equator. The 
dependence of B v on r is uncertain in the radiation zone. 
Therefore, we consider different possibilities, varying the pa- 
rameter a. The radial dependence in Eq.(8) is the simplest 
one and convenient from the computational point of view. 
A power law dependence on the cylindrical radius was also 
used, for example, in the pioneering paper by Tayler (1973) 
and, therefore, it is more convenient to compare the results 
for spherical and cylindrical geometries, choosing the depen- 
dence (8). Eq.(8) can mimic various physical situations. The 
case a > 0, for instance, can be a representative of the ra- 
diation zone in a star with the outer convective zone. On 
the contrary, the case a < can mimic a radiation zone 
with the magnetic field generated in the inner convective 
core. Certainly, the magnetic field in a radiation zone can 
have a more complicated dependence than Eq.(8) but the 
main qualitative features can be understood from this sim- 
ple model. 

It is convenient to introduce the dimensionless quanti- 
ties 



go 



77 : 



2D. 



(9) 



where ujao = B /Ry^Wp. We will represent the results in 
terms of T and 77. 

In Fig. 1, we plot the growth rate and frequency of un- 
stable modes as functions of r\ for different 9 and a = 3. 
Such distribution of B v can mimic, for example, the radia- 
tion zone of a star with a convective envelope. In this case, 
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the bottom of a convection zone is likely the location of the 
toroidal field generated by a dynamo action. However, this 
case can be also the representative of a star with relic mag- 
netic fields. Details of the formation of such fields are very 
uncertain but it is quite possible that differential rotation 
is stronger in the outer layers at the early stage of stellar 
evolution. Then, the toroidal field generated by differential 
rotation should be stronger in the outer layers as well. 

Fig. 1 shows that the instability is most efficient at the 
equator and does not occur around the pole. There always 
exist a range of 9 around the pole where the instability does 
not occur and this range depends on a. Note that the insta- 
bility exhibits a similar behaviour also in the non-rotating 
case (see Bonanno & Urpin 2012). Our result is at variance 
with a widely accepted opinion that toroidal magnetic con- 
figurations are always unstable at the axis (see, e.g., Spruit 
1999). This opinion is based on similarity of the spherical 
magnetic configuration near the axis and the axisymmetric 
cylindrical configuration. However, this analogy is generally 
incorrect because, in spherical geometry, the toroidal field 
near the axis depends also on distance along the magnetic 
axis and not only on the distance from it. This dependence 
plays the crucial role in stability. At small r/, the maximum 
growth rate is of the order of ujaq and is reached at the equa- 
tor. However, the growth rate clearly shows some suppres- 
sion for a faster rotation. Suppression becomes important 
already at relatively small values of ij ~ 1 — 2. We can distin- 
guish two substantially different regimes of suppression. If 9 
is greater than some characteristic value, do, the growth rate 
decreases with an increase of r/ approximately as 1/rj and it 
does not vanish even at very large r\. It turns out, there- 
fore, that rotation can never suppress the Tayler instability 
in the region around the equator, where 9 > 9 , but can 
only decrease its growth rate. Note that the Tayler's modes 
are complex in this region in contrast to the non-rotating 
case. The frequency (= Im T) is basically comparable to the 
growth rate and also decreases if rotation becomes faster. 
Note also that similar dependence of the growth rate on r\ 
was obtained by Spruit (1999) who considered the instabil- 
ity in cylindrical geometry with the toroidal field dependent 
on the cylindrical radius alone, — B v (s) where s is the 
cylindrical radius. The exact value of 8o is rather difficult 
to calculate because of computational problems but, in the 
case a = 3, it is equal approximately to 26°. For 9 < 9o, 
the effect of rotation is substantially stronger. In contrast to 
the instability near the equator, the instability at 9 < 6o is 
characterized by the threshold, r) cr . The instability occurs if 
?] < ri cr (Re T > in this region) but it does not occur if 
ij > rjcr since T is imaginary at such r\. The threshold, ij cr , 
depends on 9 and is ~ 1 if 9 = 24°. The threshold value 
r\ cr is lower for smaller 8. Therefore, even a relatively slow 
rotation with r\ < 1 can suppress the instability entirely in 
some region around the pole. Note that the Tayler modes 
become oscillatory beyond the threshold (at r\ > r) cr )- their 
growth rate is equal to zero but the frequency is always non- 
vanishing. In the region around the pole (at 9 < 9 cr ), the in- 
stability does not arise even if rotation is very slow (rj <C 1). 
The critical value of the polar angle, 9 cr , that restricts the 
stable region for a non-rotating star, is ~ 21°. 

In Fig. 2, we plot the growth rate and frequency for 
the case a = 2. Qualitatively, the behaviour is same but 
suppression turns out to be more essential. Again, the in- 
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Figure 1. The growth rate (left panel) and frequency (right 
panel) of the Tayler's modes as functions of the rotational pa- 
rameter 7) for a = 3 and for 6 = 90° (solid), 63° (dashed), 46° 
(dash-and-dottcd), 34° (dash-dot-dottcd), and 24° (dotted). The 
longitudinal wavenumber is / = 10. 

stability is most efficient at the equator and does not occur 
around the pole. This is qualitatively clear because, on the 
equator, with a neutral entropy gradient, we have an equiv- 
alent situation to the cylindrical geometry. The growth rate 
is maximal for small rj and is « ujao at the equator that 
is a bit lower than in the case a = 3. Rotational suppres- 
sion becomes important already at r\ ~ 1 and can operate 
in two different regimes, depending on 9. The characteris- 
tic value 9o that distinguishes these two regimes is « 37° 
in this case. Unfortunately, we cannot calculate the growth 
rate when both 9 and r\ are close to their critical values, 9o 
and ri cr , respectively. This is caused by computational prob- 
lems. At 9 > 9o, the growth rate decreases with an increase 
of r\ approximately as I/17 and it does not vanish at large 
r\. Rotation can never suppress the instability in this re- 
gion of 9 but can only reduce its growth rate. Like the case 
a = 3, unstable modes are oscillatory in this region with 
the frequency being comparable to the growth rate. Closer 
to the pole, at 9 < 9o, the instability is characterized by the 
threshold, r\ cr : the instability can occur if r\ < r\ cr but it is 
suppressed if r\ > r\ cr . For example, at 9 — 37°, the critical 
value r\ cr is w 3.5 but it is lower for smaller 9. In the region 
near the pole, 9 < 9 cr — 24°, the instability is completely 
suppressed and does not occur at any r). However, there exist 
stable oscillatory modes with a low frequency. 

Fig. 3 plots the growth rate and frequency for the 
toroidal field with a = 1. Like the previous cases, the in- 
stability is most efficient at the equator and its growth rate 
decreases if 8 decreases. Generally, the instability at a = 1 
is suppressed stronger than in the cases a = 3 and a = 2. 
However, there is a qualitative difference in comparison to 
these cases: there is no regime near the equator in which the 
growth rate is non-vanishing at large rj or, in other words, 
#o = 7r/2. It turns out that such regime of instability can 
occur only at a relatively large a. The instability at a — 1 is 
characterized by the threshold, r/ cr , even in the region near 
the equator. The threshold is lower than for a = 2 and is 
equal to « 2 at the equator. As usual, the threshold is lower 
for smaller 9. The field near the magnetic axis turns out to 
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Figure 2. The growth rate (left panel) and frequency (right 
panel) of the Tayler's modes as functions of the rotational pa- 
rameter 7) for a = 2 and for 9 = 90° (solid), 45° (dashed), 37° 
(dash-and-dottcd), 26° (dotted), and 23° (thick solid). The lon- 
gitudinal wavenumber is / = 10. 



be stable. The instability does not occur at any r\ in the 
region with 9 < 9 cr ^31°. 

The condition of stability rj > r\ cr can easily be refor- 
mulated in terms of the angular velocity and magnetic field. 
Rotation suppresses completely the Tayler instability if the 
star rotates with the angular velocity 



n > 



B 



(10) 



2 Ry/Wp' 

We can also rewrite this inequality as the condition for the 
magnetic field, 



Bo < 



(11) 



The Tayler modes become oscillatory (Im T = 0) if 
condition (10) or (11) are satisfied. Note that such behaviour 
was also seen in numerical modelling of the Tayler instability 
by Braithwaite (2006). The frequency of marginally stable 
waves is of the order of ujao{uao/Q) at r\ > 1 and it decreases 
oc at large r). This new type of magnetohydrodynamic 
waves is determined by the Coriolis and Lorentz force and 
can exist at those values of r/ which suppress the Tayler 
instability. The dispersion relation for these waves can easily 
be obtained from Eq.(6) in short wavelength approximation. 
Consider perturbations with a very short radial wavelength 
for which one can use a local approximation in the radial 
direction, such as vi r oc exp(— ik r r), where k r is the radial 
wavevector. If k r ^> max(ke, fc^), then Eq.(6) yields with 
the accuracy in terms of the lowest order in (fc r r) _1 the 
following dispersion equation 



o\ + uj\ + DQ 2 = 0, 



or- 
ai + o((2lj 2 a + Q 2 ) + cja = 0. 
The solution of this equation is 

1 



1^2 ^,2, ± l n2 /, , M 



-_(n? + a4)±-n^/i + -^ 



(12) 



(13) 



(14) 



Figure 3. The same as in Fig.l but for a = 1. The curves cor- 
respond to 9 = 90° (solid), 45° (dashed), and 37° (dash-and- 
dottcd). 



If f2; 3> uja then we can expand a square root in a power 
series of (uM/fi;) 2 - Then, choosing the upper sign in Eq.(14), 



we obtain with accuracy in the lowest order in (oM/Oi 



-W A (UA/Qi) 2 



(15) 



This dispersion relation describes a new type of oscillatory 
modes that can exist in rapidly rotating stars. These modes 
can be called the "magneto-inertial" waves because they can 
exist only in a magnetized and rotating plasma. In the con- 
sidered case, these waves are stable but, likely, they can be 
unstable under certain conditions. We consider the magneto- 
inertial waves in more detail elsewhere. 

Fig. 4 shows the growth rate and frequency for a — 
—0.4. The instability is very much suppressed by rotation in 
this case. It can occur only in a narrow region around the 
equator, 90° > 9 > 6 cr £s 65°, and does not occur in the 
extended region around the rotation axis, 9 < 6 cr « 65°. 
The threshould value of r\ at the equator is w 2. Therefore, 
the Tayler instability is entirely suppressed everywhere in 
the radiation zone with a — —0.4 if rj > 2, or fl > ujao- 
However, stable oscillating modes can exist even at much 
higher r\. 

Magnetic configurations with a rapidly decreasing 
toroidal field are stable in a cylindrical geometry. For ex- 
ample, the field can be unstable to non-axisymmetric per- 
turbations only if din B v (s)/d\ns > -1/2 (Tayler 1973a). 
Our calculations do not show the presence of instability in 
the spherical geometry if a < —1/2. 

Fig. 5 plots two examples of eigenfunctions for a = 2 
and rj — 2. The eigenfunction at 6 = 90° corresponds to 
the unstable mode that cannot be suppressed by rotation. 
Motions in this mode take place basically in the outer part 
of the radiation zone. Note that, at larger rj, this mode (cor- 
responding to 6 = 90°) turns out to be localized in a more 
narrow region near the outer boundary. On the contrary, the 
eigenfunction at 9 = 26° corresponds to a stable oscillating 
mode that exists beyond the instability threshold, rj > rj cr . 
In this case, the mode tends to be located near the inner 
boundary. This type of modes also becomes sharper with 
increasing rj. Note that this seems to be a rather general 
rule (at least, we did not find an exception): modes tend 
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Figure 4. The same as in Fig.2 but for a = —0.4 and = 90° 
(solid), and 69° (dashed). 
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Figure 5. The eigenfunctions v\ r a = 2, rj = 2, and 8 = 90° 
(solid) and 26° (dashed). 

to be localized in the outer part of the radiation zones for 
6 > 9o and, on the contrary, they are localized in the in- 
ner part for 6 < 6q. Therefore, for example, in the cases 
a = 1 or —0.4, where the region with 9 > Go is absent, both 
stable and unstable modes are localized closer to the inner 
boundary. 



4 DISCUSSION 

We have considered the stability of the toroidal field in rotat- 
ing stellar radiation zones. Stability properties of the spher- 
ical magnetic configurations are generally qualitatively dif- 
ferent from those of the cylindrical configurations. For in- 
stance, our calculations show that the instability is most 
efficient near the equator but the toroidal field can be stable 
near the symmetry axis. This is in the contrast to a widely 
accepted opinion based on the analysis of cylindrical config- 
urations (see, e.g., Spruit 1999, Zhang et al. 2003) that the 
toroidal field is always unstable near the axis in the spher- 
ical geometry. However, a direct analogy between stability 
of an infinitely long cylinder with the azimuthal field and 



the toroidal field in stellar radiation zones is generally mis- 
leading. In order to apply the cylindrical geometry near the 
axis of a spherical star, one needs to make sure that i) the 
length scale of perturbations in the axial direction is small 
compared to the radius and ii) the radial (cylindrical) length 
scale of perturbations is small compared to the radius. It is 
clear that perturbations should be affected by the spherical 
geometry if their radial length scale is comparable to the 
stellar radius. Therefore, the behaviour of such modes with 
a large radial length scale may be different in the spheri- 
cal and cylindrical geometries. As far as perturbations with 
a small radial length scale are concerned, they are alway 
stable (see Eq.(14)-(15)). Therefore, the analogy with an in- 
finitely long cylinder is generally incorrect and the Tayler 
instability can be suppressed near the magnetic axis in the 
spherical geometry. It seems that a dependence on the axial 
coordinate always makes perturbations more stable. 

Rotation provides an additional stabilizing influence on 
the Tayler instability. The effect of rotation is characterized 
by the parameter rj = 20,/loao that can be large in radiation 
zones (~ 10 2 -10 3 if B ~ 10 4 G). A reduction of the growth 
rate becomes significant already at a relatively low angular 
velocity, f2 ~ ujao, that corresponds to rj ~ 1. It turns out 
that the effect of rotation can be twofold and depends cru- 
cially on the radial dependence of the toroidal field. If the 
toroidal field increases with the spherical radius sufficiently 
rapidly (a = 2 — 3) then rotation cannot suppress completely 
the Tayler instability in the region around the equator, even 
if rj is very large. This region is more extended for larger a 
but shrinks around the equator and disappear for smaller a. 
The growth rate of instability in this region is non-vanishing 
for any angular velocity but it can be substantially reduced 
at large rj. At large rj, the growth rate descreases oc 1/rj 
and is approximately equal to u)Ao{oJAa/£l). This expression 
was also obtained by Spruit (1999) for the growth rate near 
the rotation axis in a particular case B^ = B v {s). Note, 
however, that the Alfven timescale, oJao > ^ s typically short 
(~ 1 — 3 yrs if B ~ 10 3 G) compared to the stellar life- 
time, therefore even a suppressed instability with a reduced 
growth rate can be important in radiation zones. 

Such region near the equator exists only if a is suffi- 
ciently large (> 1.5 — 2) but it shrinks and disappears if a 
decreases. At smaller 6, the region is located where the ro- 
tational suppression is qualitatively different and where the 
instability is determined by the threshold value of rj. The 
Tayler instability turns out to be suppressed completely in 
this region if rj > rj cr > 1 — 2. Therefore, modes are stable ev- 
erywhere in this region for such rapidly rotating stars. The 
threshold is not very high and corresponds approximately 
to f2 ~ u>ao- Higher eigenmodes are suppressed stronger 
than the fundamental one and perturbations with a short 
radial wavelength are always stable. Since the instability is 
suppressed at rj > rj cr but not suppressed at rj < r/ cr , this 
implies that the magnetic field should satisfy the condition 

Bo > —nR^/I^> (16) 

rjcr 

in order the instability could occur. Estimating flR ~ 2 x 10 5 
cm/s and p ~ 0.1 g/cm 3 , we obtain that instability can arise 
in the radiation zone of the Sun if Bo > 10 G This estimate 
of the critical field is more than two orders of magnitude 
higher than that obtained by Kichatinov & Riidiger (2008). 
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The hydro-dynamic motions generated by instability can Zhang K., Liao X., & Schubert G. 2003. ApJ, 585, 1124 

be important for the transport processes in radiation zones. 
Motions in the unstable modes depend generally on the po- 
lar and azimuthal wavenumbers as well as on the parameter 
j]. In real conditions, r\ is likely large (~ 10 3 — 10 4 ) and 
the radial length scale of eigenmodes is rather short. If it is 
shorter than r/l or l/m, then one can easily estimate from 
the continuity equation that the radial component of ve- 
locity is small compared to the polar and azimuthal ones. 
However, the character of motions can be different at the 
non-linear stage when the interaction between modes be- 
comes essential. 

It turns out that a new type of magneto-inertial waves 
might exist in the radiation zones of rapidly rotating stars in 
the regions where n > r/ cr and the Tayler instability is sup- 
pressed. These waves are marginally stable and their disper- 
sion relation is given by Eq.(15) at large n. The frequency of 
this oscillation decreases with an increase of Q. These waves 
are stable in our simplified model but they can be unstable 
in more realistic conditions (the presence of a poloidal field, 
differential rotation, etc.). The magneto-inertial waves can 
play an important in various processes in the radiation zone 
(mixing, transport of the angular momentum, formation of 
the tachocline, etc.). These waves will be considered in more 
detail in a subsequent publication. 
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